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ABSTRACT 

The spectral characteristics of the pulsed gamma-ray emission from outer-magnetospheric particle 
accelerators are investigated. Either positrons or electrons are accelerated outwards by the magnetic- 
field- aligned electric field to emit gamma-rays via curvature process. Since the particles move along 
relatively straight paths in the trailing side of a rotating magnetosphere, they attain higher Lorentz 
factors to emit more energetic gamma-rays than those in the leading side. It is first demonstrated that 
the cutoff energy of the curvature radiation evolves with the rotation phase owing to the variation of 
the curvature radii of the particle paths and maximizes at a slightly later phase of the trailing peak 
in the gamma-ray light curve. 

Subject headings: gamma rays: stars — magnetic fields — methods: analytical — methods: numerical 
— stars: neutron 



1. INTRODUCTION 

The Fermi Large Area Telescope (LAT) provides a 
wealth of new data on isolated, rotation-powered pul- 
sars, increasing the number of detected 7-ray pulsars 
from seven to more than sixty (e.g., Abdo et al. 2010a; 
Abdo et al. 2010b; Saz Parkinson et al. 2010). The 
Astro-rivelatore Gamma a Immagini LEggero (AGILE) 
has also reported the detection of about twenty 7-ray pul- 
sars (e.g., Pellizzoni et al. 2009; Pilia et al. 2010). Since 
interpreting 7-rays is less ambiguous compared with re- 
processed, non-thermal X-rays, the 7-ray pulsations ob- 
served from these objects are particularly important as a 
direct signature of basic non-thermal processes in pulsar 
magnctosphcres, and should help to discriminate among 
different emission models. In this letter, we thus focus on 
the spectral evolution of the pulsed 7-rays and examine 
if the hardening of the trailing peak in the light curve, 
a general property derived from these observations, can 
be explained in the context of the outer-magnetospheric 
emission scenario. 

The light curves and spectral evidence obtained by 
the observations mentioned above, suggest that the 7- 
ray pulsars have high-altitude emission zones whose fan- 
like beams scan over a large fraction of the celestial 
sphere. Therefore, recent pulsar high-energy emission 
models adopt higher-altitude emission geometries. There 
are three main scenarios in this approach: the outer-gap 
(OG) model (Cheng et al. 1986a, b; Romani 1996; Cheng 
et al. 2000; Hirotani 2006b, 2008; Takata et al. 2008; 
Romani & Watters 2010), the higher-altitude slot-gap 
model (Muslimov & Harding 2004), and the pair-starved 
polar-cap model (Frackowiak & Rudak 2005; Harding et 
al. 2005; Venter et al. 2009). 

In the present letter, we focus on the OG model, 
in which pairs are created in the outer magnetosphere 
mainly by photon-photon (7-7) pair production. The 
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produced pairs polarize owing to the magnetic-field- 
aligned electric field, En . If the rotation and mag- 
netic axes reside in the same hemisphere, a positive -Ey 
is exerted to accelerate e + 's (or e~'s) outward (or in- 
wards), increasing the real charge density outwards. The 
inward-migrating, relativistic e~'s radiate curvature 7- 
rays, some of which (nearly head-on) collide with the X- 
rays emitted from the NS surface to materialize as pairs. 

In a pulsar magnetosphere, there is a surface called 
the 'null surface', on which the Goldreich- Julian charge 
density (Goldreich & Julian 1969) p G3 = —B ■ n/(2irc) 
changes sign, where B denotes the local magnetic field, f2 
the NS angular-velocity vector, and c the speed of light. 
There is another characteristic surface called the 'light 
cylinder' beyond which the co-rotational velocity exceeds 
c. The distance of the light cylinder from the rotation 
axis is called the 'light-cylinder radius', m^c = c/fi, 
where = \Sl\ denotes the NS rotational angular fre- 
quency. On each magnetic azimuthal angle (measured 
around the magnetic axis counter-clockwise), there is a 
magnetic field line that crosses the light cylinder tangen- 
tially; they are called the 'the last-open field lines'. An 
OG is essentially located between the null surface and 
the light cylinder in the higher colatitudes (i.e., in the 
closer region to the magnetic axis) than the last-open 
field lines. 



2. CHARACTERISTIC ENERGY OF CURVATURE 
RADIATION 

The characteristic energy of curvature emission is given 
by (e.g., Rybicki & Lightman 1979) 
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where 7 refers to the Lorentz factor of the electron or 
positron, h the Planck constant divided by 2w, and g c 
the curvature radius of particle's three-dimensional (3- 
D) motion. In the gap, the e ± 's achieve electrostatic 
force balance, eEu = 2e 2 j i /(3gl), and saturate at the 
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terminal Lorentz factor, 
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where e designates the charge on the positron, and Eu 
the electric field component projected along the magnetic 
field line. Combining equations (JIJ and ©, we obtain 
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Since 



are saturated at 7, the spectral den- 
sity of their curvature emission declines sharply as 
{E/E c ) 1/2 exp(-E/E c ) above the cutoff energy, E > E c . 
Therefore, for the same En , equation ([3]) shows that the 
curvature spectrum has a greater E c for a greater g c . 

Since the pair production mainly takes place between 
the inward 7-rays and the surface X-rays, most of the 
pairs are produced in the lower altitudes where the X- 
ray density is large and the collisions take place mostly 
head-on. Thus, outward- moving species (e.g., e + 's) run 
longer distances than inward ones within the gap, lead- 
ing to a stronger outward 7-ray flux than the inward one. 
Therefore, to elucidate the physical mechanism that pro- 
duces harder spectrum in the trailing peak of a 7-ray 
light curve, we concentrate on the outward-moving par- 
ticles and the resultant outward emissions in § [5] 

To reveal the nature of curvature radiation, we must 
consider the curvature radius of actual particle motion 
in the 3-D pulsar magnetosphere. Well inside the light 
cylinder, the drift motion due to magnetic curvature, 
magnetic gradient, or time- varying electric field is small 
compared to the E x B drift. Thus, particles co-rotate 
with the magnetic field lines. In the flaring field line 
geometry of a NS magnetosphere, the path of outward- 
moving particles can be illustrated as figure [1] In the 
leading side, particle's path is bent towards the rotational 
direction as the thin solid curve indicates. However, in 
the trailing side, the path is straighter as indicated by the 
thick solid curve. Thus, outward-moving particles have 
greater curvature radii in the trailing side than in the 
leading side. This leads to a greater E c in equation 
This is the qualitative reason why most of the 7-ray pul- 
sars detected with LAT exhibits harder spectrum in the 
trailing light-curve peak than in the leading peak. 

Next, let us examine the distribution of the curvature 
radius more quantitatively. Due to strong relativistic 
beaming, curvature photons are emitted along the in- 
stantaneous particle velocity measured by the distant 
static observer (i.e., us). In the polar coordinates, the 
orthonormal components of the particle velocity become 
(Mestel et al. 1985; Camenzind 1986a,b) 
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J {B r ) 2 + (B 6 ) 2 , B denotes the strength of the 
magnetic field, m the distance from the rotation axis. 




Fig. 1. — Schematic figure (top view) of particle motion in a ro- 
tating pulsar magnetosphere. The thick dashed lines represent the 
dipolar magnetic field lines at some moment, and the thin dashed 
ones do the same field lines at later rotation phase. The thin and 
thick solid lines represent the paths of particles moving outwards 
in the leading side (LS) and trailing side (TS), respectively. 



Both {v r ,v e ,v*) and (B r ,B 6 ,B&) are measured by a dis- 
tant static observer (i.e., in the inertial observers frame). 
The upper sign of f v (eq. [3]) corresponds to the outward 
motion along the magnetic field line, whereas the lower 
sign to the inward motion. These expressions are quite 
general and valid irrespective of the force balance on the 
particles as long as 7 3> 1, because they are derived only 
from the Maxwell equations and the frozen-in condition, 
the latter of which is justified since we consider En <C B 
(or equivalently, consider that the potential drop in the 
gap is much small compared to the electro-motive force 
exerted on the spinning NS surface). For example, they 
are valid not only for e s being accelerated without dis- 
sipation but also for those in the force balance between 
radiation-reaction and electro-static acceleration. Note 
that equations (j3| and ([5]) correspond to a generalization 
of § 4.1 of Dyks et al. (2010) into the higher altitudes. 
Note also that equations Q and ([5]) are applicable not 
only within the light cylinder but also outside of it and 
that < holds in ordinary situations. 

Integrating equation (j4|) along particle motion, we ob- 
tain the path of the particle, which allows us to compute 
q c as a function of the distance s along the field line. 
Since the incorporation of a field-line deformation due to 
magnetospheric currents is out of the scope of this letter, 
we adopt a retarded, vacuum dipole field (Cheng et al. 
2000) as an example. In figure [21 we present q c /wlc as 
a function of s for a representative magnetic field line in 
the leading side (LS; = +60°) and that in the trail- 
ing side (TS; (p* = —60°), where tp* refers to the mag- 
netic azimuthal angle measured around the magnetic axis 
counter-clockwise. The magnetic inclination angle, a, is 
assumed to be 60°. 

In figure [3J the (red) thick solid and dashed curves 
denote Qc/^lc of the actual particle path, whereas the 
(black) thin solid and dashed curves do £ c /^lc of the 
last-open magnetic field lines. The LS quantities are plot- 
ted with dashed curves, while the TS quantities with 
solid curves. It follows that g c /^LC of the particle path 
(labeled with P-TS) becomes greater than that of the 
field line (B-TS) in the TS. In the LS, on the contrary, 
£>c/wlc the particle path (P-LS) becomes less than 
that of the field line (B-LS). Note that B-TS is corrected 
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Fig. 2. — Distribution of the curvature radius, q c as a function 
of the distance along the magnetic field line. Both the abscissa 
and the ordinate are normalized by the light cylinder radius, tolc- 
The (red) thick solid curve (labeled as P-TS) denotes £> c /roLC °f 
the actual particle path in the TS (specifically, at ip„ = —60°). 
The (red) thick dashed curve (P-LS) denotes Qc/i&lc of the actual 
particle path in the LS (<p* = 60°). The (black) thin solid curve 
(B-TS) does gc/^LC °f the last-open magnetic field line in the TS 
(</?* = —60°), while the (black) thin dashed curve (B-LS) that of 
the last-open magnetic field line in the LS (yj* = 60°). 

to P-TS and B-LS to P-LS, and that the magnetic field 
lines (see B-LS vs. B-TS) have greater g c in the LS 
whereas the particle paths (P-LS vs. P-TS) have greater 
g c in the TS. In previous pulsar emission models, it has 
been phenomenologically assumed that the particles emit 
curvature radiation along local magnetic field lines and 
adopted B-LS and B-TS as g c /m-^c', thus, it has been im- 
possible to reproduce a harder spectrum in the TS than 
the LS. Precisely speaking, they did considered the aber- 
ration of photon momentum to compute the emission di- 
rection from each point (and to reproduce the observed 
pulse profiles) , but have computed g c from the magnetic 
field geometry (as represented by the thin lines labeled 
with B-LS and B-TS). Instead, we should take account 
of particle's' motion with equation ((4]) and adopt the 
correct g c (as represented by P-LS and P-TS); then, we 
obtain a harder spectrum in the TS as observed from 
most of the 7-ray pulsars. 

3. SELF-CONSISTENT ACCELERATOR SOLUTION 

To precisely compute the phase-resolved spectrum, we 
must solve En and the distribution functions of e s at 
each point in the 3-D pulsar magnetosphere, together 
with the radiative transfer. In this section, we thus solve 
the set of Maxwell and Boltzmann equations in the 3- 
D pulsar magnetosphere in the open zone from the NS 
surface to the light cylinder under appropriate bound- 
ary conditions (see Hirotani 2006a, b for details). To 
demonstrate general spectral properties, we consider typ- 
ical young NS parameters, P — 0.1 s, P = 10~ 13 s s" 1 , 
kT s = 50 eV, where kT refers to the blackbody tempera- 
ture of the whole NS thermal emission. For this 'example 
pulsar', we adopt a — 60°. 

We present the distribution of the gap in the magne- 
tosphere by showing the gap fractional thickness, /i m , 
projected on the last-open field line surface in figure [3J 
If h m = 1.0, it means that all the open fluxes threads 
the OG. On the other hand, if h m <C 1, only the fluxes 
very close to the last-open field lines are active. It fol- 
lows that the gap exist mainly outside the null surface 
(white curve); thus, the present solution coincidentally 
corresponds to a quantitative extension of previous, phe- 



nomenological OG models. 

We next present the distribution of photon emission 
as a function of NS rotational phase and the observers' 
viewing angle £ in figure01 where the distance is assumed 
to be 1 kpc. Since only the photons emitted from the OG 
connected to the north magnetic pole are plotted, pho- 
tons propagating into the southern hemisphere (( > 90°) 
are mostly emitted by outward-moving charges, whereas 
photons into £ < 90° are by inward ones. Thus, we find 
that the outward photon flux dominates the inward one. 

We can find the pulse profiles by slicing across the pho- 
ton intensity plot in figure U at a constant f. As an ex- 
ample, we choose £ = 115° and show the pulse profile in 
figure [5] as the dashed line (in arbitrary unit); PI corre- 
sponds to the leading peak, while P2 the trailing peak. 
One full NS rotation period is depicted in the abscissa. 

Finally and most importantly, we examine the 
evolution of spectral properties with phase. To 
this aim, we fit the solved differential photon flux 
(photons cm~ 2 s _1 MeV -1 ) with a power law with ex- 
ponential cutoff shape given by 
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in each phase bin with 1/72 (i.e., 5°) interval, where both 
E and E c are in GeV unit and refer to the photon energy 
and the cutoff energy, respectively; T is the photon index. 
The best fit T and E c are depicted in figure [5j The (red) 
thick solid curve in the top panel denotes the evolution 
of r with phase, while that in the bottom panel that of 
E c . For comparison, T and E c that would be obtained 
if we adopted p c of the magnetic field lines (instead of 
that of the particle path) are also depicted by thin dotted 
curves. It follows that the photon index remains roughly 
constant during the bright phase, while the cutoff energy 
sharply peaks in the trailing peak as expected from the 
analytic argument in § [2j Note that only comparable 
E c would be obtained in the two light-curve peaks if we 
adopted p c of the magnetic field lines (as in previous 
pulsar emission models) . The conclusion of greater E c in 
the trailing peak is, in deed, consistent with the spectral 
properties of most of the 7-ray pulsars detected with LAT 
(for PSR J2021+3651, see Abdo et al. 2009; for the Crab 
pulsar, see Abdo et al. 2010b; for the Vela pulsar, see 
Abdo et al. 2010c; for PSR J1057-5226, PSR J1709- 
4429, PSR J1952+3252, see Abdo et al. 2010d; for the 
Gcminga pulsar, see Abdo et al. 2010e). 

4. DISCUSSION 

By examining the actual paths of outward-moving 
charged particles in a 3-D pulsar magnetosphere, we 
demonstrated that the curvature radii of their paths be- 
come larger in the trailing side (TS) than in the leading 
side (LS) in § Applying this result to the curvature 
radiation formula (§ [l}, we first gave a physical expla- 
nation why the trailing peak generally exhibits a harder 
spectrum than the leading peak in 7-ray light curves. 
To confirm this analytical prediction, we solved the set 
of Maxwell and Boltzmann equations in the 3-D pulsar 
magnetosphere in § [31 finding that the obtained emis- 
sivity distribution does result in a greater cutoff energy 
in the trailing peak than in the leading peak. Although 
the discussion in § [3] is limited to the OG model, the 
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Fig. 3. — Distribution of the gap meridional thickness projected 
on the last-open field line surface. The white curve denotes the 
cross section between the last-open field lines and the null surface. 
The ordinate denotes the angle, [rad], measured around the 
magnetic axis counter-clockwise. The region ip* > (or ip„ < 0) is 
referred to as the leading (or trailing) side. 
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Fig. 4. — Angular distribution of photon emission above 89 MeV 
(log 10 of flux per phase in MeVs~ 1 cm -2 deg -1 unit) from the 
outer gap connected to the north magnetic pole. The phase 0° 
(in the abscissa) corresponds to the neutron-star rotation phase 
at which the photon emitted from the north polar cap arrives the 
observer. The ordinate denotes the observer's viewing angle f with 
respect to the rotation axis. 

same argument can be, in fact, readily applied to other 
higher-altitude emission models, as long as the emission 
is due to the curvature process. What we have to change 
is to amend E\\ (instead of what is obtained in the OG 
model) as a function of position in the 3-D pulsar magne- 
tosphere. Unless En becomes substantially greater in the 
LS than in the TS, the same conclusion will be derived, 
which is, in deed, expected both in the slot-gap and in 
the pair-starved polar-cap models. 

In figure [3 T peaks at pulsar phase ~ —30°. However, 
this peaking structure of T appears at small-flux phase, 
and hence is vulnerable to the choice of £. For example, 
it disappears for £ > 120°. On the contrary, the peaking 
structure of E c at a slightly later phase of the trailing 
peak, is generally obtained in a wide parameter range 
of (P,P,kT,a) and £. Although not all the parameter 
space has not been investigated, the conclusion that the 
trailing peak has a greater cutoff energy than in other 
phases, appears to be universal. 

In this letter, we have adopted a retarded, vacuum 
dipole magnetic field. However, in general, we should 
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Fig. 5. — Evolution of the photon index T (top panel) and the 
cutoff energy E c (bottom panel) with the rotational phase, is de- 
picted by the (red) thick solid curve, for an example pulsar with 
P = 0.1 s, P = 10~ 13 s s~\ kT = 50 eV, and a = 60°. Observer's 
viewing angle is assumed to be f = 115°. For comparison, T and 
E c that would be obtained if we adopted the curvature radii of 
the magnetic field lines (instead of the particle paths) are depicted 
by the thin dotted curves. Pulse profile is overlaid as the dashed 
curves (in arbitrary unit). 

take account of the deformation of the field lines due to 
the magnetosphcric current (Contopoulos et al. 1999; 
Gruzinov 2005; Spitkovsky 2006; Bai & Spitkovsky 
2010a, b; Contopoulos & Kalapotharakos 2010). The 
effect of deformation works to straighten the field lines 
in the outer magnetospheres, pushing the OG position 
towards the light cylinder particularly in the TS; for ex- 
ample, if a = 60°, the OG within the region (/j* < —60° 
is suggested to be less active due to the position shift to- 
wards the light cylinder. Nevertheless, the self-consistent 
OG solution presented in §[3] shows that the pair produc- 
tion is not, any way, self-sustained in ip* < —60°, result- 
ing in a negligible emission from this region. For exam- 
ple, the trailing peak (around —10° in pulsar phase) is 
mostly formed by the emission from —60° < tp* < —30°. 
In addition, the specific intensity towards £ ~ 115° di- 
rection will not be strongly affected by the straighten- 
ing, which follows from panel (d) of figure 8 in Bai and 
Spitkovsky (2010b). On these grounds, we consider that 
the conclusions in the present letter are relatively less 
vulnerable to the correction of the field-line geometry 
due to magnetospheric currents. 

We can further apply the present numerical method to 
any rotation-powered pulsars. In subsequent papers, we 
will demonstrate that the numerical solutions are consis- 
tent with various observations of individual pulsars, such 
as the energy fluxes, pulse profiles, and phase-resolved 
spectra. To this aim, it is essential to self-consistently 
solve the screening processes due to the polarization of 
the produced pairs, together with the distribution func- 
tions of e s and photons. This can be achieved only by 



Hardening of Trailing Peak in Gamma-ray Pulsar Light Curve 



5 



solving the Maxwell and Boltzmann equations simulta- 
neously in the 3-D pulsar magnetosphere. 
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